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Abstract. The possibility of determining cosmological 
parameters on the basis of a wide set of observational 
data including the Abell-ACO cluster power spectrum and 
mass function, peculiar velocities of galaxies, the distribu- 
tion of Ly-a clouds and CMB temperature fluctuations 
is analyzed. Using a y 2 minimization method, assuming 
J7a + f2 m atter = 1 and no contribution from gravity waves, 
we show that this data set determines quite precisely the 
values of the spectral index n of the primordial power 
spectrum, baryon, cold dark matter and massive neutrino 
density f2(,, Q c dm and f2„ respectively, the Hubble constant 
h = i/o/(100km/s/Mpc) and the value of the cosmological 
constant, Q\ • 

Varying all parameters, we found that a tilted AMDM 
model with one sort of massive neutrinos and the parame- 
ters n = 1.12±0.10, n m = 0.41 ±0.11 (0 A = 0.59±0.11), 
n cdm = 0.31±0.15, Q„ = 0.059±0.028, fl b = 0.039±0.014 
and h = 0.70 ± 0.12 matches observational data best. 

fli, is higher for more species of massive neutrinos, ~ 
0.1 for two and ~ 0.13 for three species. ft m raises by 
~ 0.08 and ~ 0.15 respectively. 

The la (68.3%) confidence limits on each cosmological 
parameter, which are obtained by marginalizing over the 
other parameters, are 0.82 < n < 1.39, 0.19 < fi m < 1 
(0 < Q A < 0.81), < n v < 0.17, 0.021 <n b < 0.13 and 
0.38 < h < 0.85 1.5 < b c i < 3.5. Here b c i is the cluster bias 
parameter. The best-fit parameters for 31 models which 
are inside of la range of the best model are presented 
(Table 4). 

Varying only a subset of parameters and fixing the 
others changes the results. In particular, if a pure matter 
model (fl m = 1) is assumed, MDM with fl v = 0.22 ±0.08, 
three species of massive neutrinos and low h = 0.47 ±0.05 
matches the observational data best. If a low density 
Universe Q m = 0.3 is assumed, a ACDM model with- 
out hot dark matter and high h = 0.71 matches the 
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observational data best. If the primordial power spec- 
trum is scale invariant (n = 1) a low density Universe 
(fi m = 0.45 ± 0.12, h = 0.71 ± 0.13) with very little hot 
dark matter (f2„ = 0.04 ± 0.03, N v = 1) becomes the best 
fit. 

It is shown also that observational data set used here 
rules out the class of CDM models with h > 0.5, scale in- 
variant primordial power spectrum, zero cosmological con- 
stant and spatial curvature at very high confidence level, 
> 99.99%. The corresponding class of MDM models are 
ruled out at ~ 95% C.L. 



Key words: Large Scale Structure: cosmological models, 
power spectrum, cosmological parameters 

1. Introduction 

Observations of the large scale structure (LSS) of the Uni- 
verse carried out during the last years and coming up from 
current experiments and observational programs allow to 
determine the parameters of cosmological models and the 
nature of dark matter more precisely. The usual cosmolog- 
ical paradigm - a scale free power spectrum of scalar pri- 
mordial perturbations which evolve in a multicomponcnt 
medium to form the large scale structure of the Universe 
- is compatible with the observed cosmic microwave back- 
ground (CMB) temperature fluctuations. Most inflation- 
ary scenarios predict a scale free primordial power spectra 
of scalar density fluctuations P{k) ~ k n with arbitrary n 
as well as gravity waves which contribute to the power 
spectrum of CMB temperature fluctuations (^j(-)t at low 
spherical harmonics. But models with a minimal number 
of free parameters, such as the scale invariant (n = 1) 
standard cold dark matter model (SCDM) or the stan- 
dard mixed (cold plus hot) dark matter model (SMDM) 
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only marginally match observational data. Better agree- 
ment between predictions and observational data can be 
achieved in models with a larger number of parameters: 
cold dark matter (CDM) or mixed dark matter (MDM) 
with baryons, a tilted primordial power spectra, spatial 
curvature (flfc), a cosmological constant (CIa) and a ten- 
sor contribution to the CMB anisotropy power spectrum. 

The neutrino oscillations discovered recently in the 
Super-Kamiokande experiment (Fukuda et al. 1998) show 
that at least one species of weakly interacting neutri- 
nos have non-zero rest mass. Assuming that the larger 
one of them is about m v ~ '5m% w 0.07 eV we find 
VL V w 7.4 x W~ 4 N v /h 2 . It is also possible that one, two 
or three species have masses in the eV range and give ap- 
preciable contribution to the dark matter content of the 
Universe. 

The presence of rich clusters of galaxies at z « 
0.54, 0.55, 0.8 (Bahcall & Fan 1998) indicates a low mat- 
ter density. 

In this work we do not include into our anal- 
ysis the recent observations of distant supernovae 
(Perlmutter et al. 1998, Riess et al. 1998). The SNela 
measurements support a positive cosmological con- 
stant. Assuming a flat Universe, Oa + = 1, a 
value of fl\ ~ 0.7 is preferred (see also the re- 
view by Bahcall et al. 1999), but, in agreement with 
Valdarnini et al. 1998, Primack & Gross 1998, we find 
that on the basis of LSS data alone, a non-vanishing cos- 
mological constant is preferred within the class of models 
analyzed in this work. 

Another approach based on the search of best-fit cos- 
mological parameters in open and critical density CDM 
and ACDM models without gravitational waves for the to- 
tal combination of observational data on CMB anisotropy 
has been carried out by Lineweaver & Barbosa 1998. But 
the CMB data set corresponds to very large scales 
(> 100/i _1 Mpc) and it is not sufficiently sensitive to 
the existence of a HDM component. The power spec- 
tra of density fluctuations obtained from the spatial 
distribution of Abell-ACO clusters (Einasto et al. 1997, 
Retzlaff et al. 1998), APM, CfA and IRAS galaxy surveys 
(Einasto et al. 1999 and references therein) are extended 
to smaller scales up to galaxy scales which are below the 
neutrino free streaming scale. On small scales constraints 
are obtained from absorption features in quasar spectra 
known as the Ly-a forest (Gnedin 1998, Croft et al. 1998). 

The determination 
of cosmological parameters from some observations of the 
LSS of the Universe was carried out in many papers (e.g. 
Atrio-Barandela et al. 1997, Lineweaver & Barbosa 1998, 
Tegmark 1999, Bridle et al. 1999, Novosyadlyj 1999 and 
references therein). Recently, Bridle et al. 1999 have an- 
alyzed the cluster abundances, CMB anisotropies and 
IRAS observations to optimize the four parameters 
(0 TO , h, as, and &iras m a open CDM model. 
Atrio-Barandela et al. 1997 use the cluster power spec- 



trum together with data of the Saskatoon experiment to 
discuss the possible existence of a built-in scale in the 
primordial power spectrum. In this paper a total of 23 
measurements from sub-galaxy scales (Ly-a clouds) over 
cluster scales up to the horizon scale (CMB quadrupole) 
are used to determine seven cosmological parameters. 

Clearly, it is possible that the 'correct cosmological 
model' is not one of those analyzed in this paper. If the 
data are good enough this can in principle be decided by 
a x 2 -test. As long as we find a model within the family 
of models studied here with an acceptable value of x 2 , we 
have no compelling reason to consider other models. 

In view of the growing body of observational data, we 
want to discuss the quantitative differences between the- 
ory and observations for the entire class of available mod- 
els by varying all the input parameters such as the tilt of 
the primordial spectrum, n, the density of cold dark mat- 
ter, fl c dm, hot dark matter, Q„, and baryons, the vac- 
uum energy or cosmological constant, Q\, and the Hubble 
parameter h, to find the values which agree best with ob- 
servations of LSS on all scales (or even to exclude a whole 
family of models). Here we restrict ourselves to the analy- 
sis of spatially flat cosmological models with + Q m = 1 
(fife = 0), where TO = Q c dm + ^6 + and to an in- 
flationary scenario without tensor mode. We also neglect 
the effect of a possible early reionization which could re- 
duce the amplitude of the first acoustic peak in the CMB 
anisotropy spectrum. 

The reason for the restriction of flat models is mainly 
numerical. However, the new CMB anisotropy data from 
the Boomerang experiment actually strongly favors spa- 
tially flat universes (Melchiorri et al. 1999b). Neglecting 
the tensor mode which affects the normalization and the 
height of the first acoustic peak is motivated by the work of 
Tegmark 1999, who found that CMB anisotropy data pre- 
fer no or a small tensor component, however there are also 
arguments in favor of the importance of the tensor mode 
(Arkhipova et al. 1998, Melchiorri et al. 1999a). Further- 
more, since the LSS data used in this paper disfavors very 
blue spectra, the high acoustic peak indicates that reion- 
ization cannot be substantial for the class of models ana- 
lyzed in this paper. Hence we set the optical depth r = 0. 

The outline of this paper is as follows: In Sect. 2 we de- 
scribe the observational data which are used. The method 
of parameter determination and some tests are described 
in Sect. 3. We present the results obtained under differ- 
ent assumptions about the parameter ranges in Sect. 4. A 
discussion of our results and the conclusions are given in 
Sects. 5 and 6 respectively. 

2. The experimental data set 

2.1. The Abell-ACO cluster power spectrum 

One might expect that the most favorable data for the 
determination of cosmological parameters are power spec- 
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tra constructed from the observed distribution of galax- 
ies. But the power spectra of galaxies obtained from the 
two-dimensional APM survey (e.g. Maddox et al. 1996, 
Tadros & Estathiou 1996, and references therein), the 
CfA redshift sur- 

vey (Vogeley et al. 1992, Park et al. 1994), the IRAS sur- 
vey (Saunders et al. 1992) and/or from the Las Campanas 
Redshift Survey (da Costa et al. 1994, Landy et al. 1996) 
differ both in the amplitude and in the behavior near 
the maximum. Moreover, nonlinear effects on small 
scale must be taken into account in their analysis. For 
these reasons we do not include galaxy power spec- 
tra for the determination of parameters in this work. 
Here, we use the power spectrum of Abell-ACO clusters 
(Einasto et al. 1997, Retzlaff et al. 1998) as observational 
input. This power spectrum is measured in the range 
0.03/i/Mpc< k < 0.2/i/Mpc. The cluster power spectrum 
is biased with respect to the dark matter distribution. We 
assume that the bias is linear and scale independent in the 
range of scales considered. The position of the maximum 
{k m ax ~ 0.05/i/Mpc) and the slope at lower and larger 
scales are sensitive to the baryon content f^, the Hub- 
ble constant h, the neutrino mass m v and the number of 
species of massive neutrinos N v (Novosyadlyj 1999). The 
Abell-ACO cluster power spectrum PA+ACo(kj) (here 
and in the following a tilde denotes observed quantities) 
has been taken from Retzlaff et al. 1998. We present 13 
values of pA+Aco(kj) and the la errors in Table 1 and 
in Fig. 1. In a first step, we have assumed that the 13 
points in this power spectrum given below are indepen- 
dent measurements. The value of x 2 obtained under this 
assumption is much smaller than the number of degrees of 
freedom (see below). We interpret this as a hint that the 
13 points of Pa+aco given in Table I cannot be consid- 
ered as independent measurements. We therefore describe 
the power spectrum by three parameters A, fcbcnd and a 
to be of the form 



PA+ACo(k) 



Ak 



1 + (fc/fcbcnd)" ' 



(1) 



A fit of the parameters to the observed power spectrum 
gives 

A = (3.78 ± 1.71) x 10 6 , fc bcnd = 0.056 ± 0.015, 



a = 3.49 ±0.72. 



(2) 



In Fig. 1 we show the observed power spectrum together 
with the fit. The cosmological model parameters obtained 
using the full power spectrum information or the three 
parameter fit are in good agreement, but the latter pre- 
scription leads to a more reasonable value of x 2 - 

This point is quite important since it illustrates that 
a small \ 2 need not mean that the error bars of the data 
are too large but it can be due to data points depending 
only on a few parameters and therefore not being inde- 
pendent. If a power spectrum, like the one above can be 



modeled by 3 parameters, then, by varying three cosmo- 
logical parameters, like e.g. the cluster bias b c , the HDM 
contribution f2„ and the Hubble parameter, h, we can in 
general (if there is no degeneracy) fit all three parameters 
A, fcbend and ol and thereby the entire power spectrum. 
The number of degrees of freedom in such a fit is and 
not 10 as one would infer form the number points of the 
power spectrum. 

To make best use of the observational information, we 
nevertheless use the full 13 points of the power spectrum 
to fit the data, but we assign it rip = 3 for the number of 
degrees of freedom. 
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Fig. 1. The Abell-ACO power spectrum by Retzlaff et al. 1998. 
The solid line is the best fit according to Eqs. 1 and 2. 
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2.2. CMB data 

We normalize the power spectrum using the COBE 4-year 
data of CMB tempera- 

ture fluctuations (Bennett et al. 1996, Liddle et al. 1996, 
Bunn and White 1997). We believe that using all avail- 
able experimental data on AT/T on angular scales smaller 
than the COBE measurement is not an optimal way for 
searching of best-fit parameters because some data points 
in CMB spectrum contradict each other. Therefore, we use 
only the position and amplitude of the first acoustic peak 
derived from observational data as integral characteristics 
of CMB power spectrum, which are sensitive to some of 
the model parameters. 

To determine the position, £ p , and amplitude, A p , of 
the first acoustic peak we use the set of observational 
data on CMB temperature anisotropy given in Table 2 
(altogether 51 observational points). For each experiment 
we include the effective harmonic, the amplitude of the 
temperature fluctuation at this harmonic, the upper and 
lower error in the temperature, and the effective range of 
the window in l-spa.ee. In those cases when original pa- 
pers do not contain effective harmonics and band width 
we have taken them from Max Tegmark's CMB data 
analysis center (Tegmark 1999) dated Nov 25 1999. We 
fit the experimental data points by a polynomial of 6- 
th order using the Levenberg-Marquardt method to de- 
termine the position and amplitude of the first peak: 
[1(1 + l)C i /27r 2 ] 1 / 2 = £® =0 a,l\ The best-fit values of the 
coefficients are: ao = 31.1, a\ = —0.309, a 2 = 5.18 x 10~ 3 , 
a 3 = 9.66 x 10~ 6 , a 4 = -3.68 x 10~ 8 , a 5 = 1.22 x 10~ 10 , 
a 6 = -9.19 x 10~ 14 (x 2 = 62.9). The amplitude A p and 
position l p of first acoustic peak determined from data 
fitting curve are 79.6^iK and 253 correspondingly. Our re- 
sult differs only slightly from the numbers obtained by 
Lineweaver & Barbosa 1998 who found 260 and 88/xK. 
Fig. 2 shows the observational data used together with 
the polynomial best fit (solid line). 

We estimated the error of A p and l p in the following 
way. 

By varying of all coefficients at we determine x 2 - 
hypcr-surfacc in 7-dimension parameter space which con- 
tains deviations of less than la. If the probability dis- 
tribution obeys Gaussian statistics, this corresponds to a 
68.3% confidence level. It is well known, that present CMB 
anisotropy data even on small scales do not obey Gaussian 
statistics and thus this procedure is somewhat arbitrary. 
However, it can be assumed that this gives us a good in- 
dication for the errors bar in position and amplitude of 
the first acoustic peak. For 44 degrees of freedom (51 data 
points minus 7 parameters) this hyper-surface corresponds 
to Ax 2 = 47.9. For values of parameters at which have a 
X 2 < Xmin + ^x 2 we calculate the Ci 's and find the peak 
amplitude A p and the position l p . They are in the contour 
line in the A p — l p plane shown in Fig. 3. The upper-lower 
and right-left extremal points indicate la statistical er- 



rors: AAf = +17.0, -16.3/jK and Al s p l = +28, -22. Un- 
certainties of effective harmonics of each experiment do 
not influence the error of the amplitude of the first acoustic 
peak but must be take into account for the full error in the 
peak position, so that Al p = Al p l + Al p , where last term 
is the mean band width around l p . We estimate it as mean 
width of all experiments weighted by acoustic peak ampli- 
tude Al p = Xw=i(A0i w i/Si=i Wj, where the weighting 
factor LVi = + 1)CiJ2tt 2 ] 1 / 2 /A p is calculated using 
polynomial fit. This finally leads to Al p w 45.0. (Without 
weighting the value is Al™ w 42). Therefore, the errors of 
determination of first acoustic peak amplitude and posi- 
tion are AA p w 16.5£iK and Al p w 70 respectively. We 
use these errors below in our search of cosmological pa- 
rameters. 

It is interesting to note that no 6th order polynomial 
fits the data really well. For our best fit polynomial we 
obtain \ 2 ~ 62.9 for 51 data points and 7 parameters. 
The probability for this polynomial leading to the ob- 
served data is about 1%. This big x 2 can have two ori- 
gins. First, the probability distribution is non-Gaussian 
and, therefore, the probability to obtain this value of x 2 
is higher than 1%. Secondly, some data seems to be con- 
tradictory. For example, if we ignore all the Python V 
points we obtain a best fit polynomial with x 2 = 23 which 
is even slightly too low. (Removing of these points does 
not change essentially the result amplitude and position 
of acoustic peak, l p = 256, A p = 79.0 without them). But 
of course we are not allowed without any good reason, 
to leave away some experimental results. It may well be 
that Python V is correct and some other experiments are 
wrong. Therefore, we adopted this somewhat hand wav- 
ing way to extract information from this data. Clearly, 
a more thorough analysis with true, non-Gaussian likeli- 
hood functions would be in order, which we leave for the 
future (see Bartlett et al. 1999). 

For the comparison of models with the CMB data we 
use, apart from the COBE normalization, only the two 
parameters obtained by the fitting procedure described 
above: the effective harmonic l p = 253 + 70 of the peak po- 
sition and the amplitude of the peak A p = 79.6 ± 16.5^K. 
Clearly, this position and height of the first acoustic peak 
is not strictly implied by the present data and can there- 
fore be criticized. In this sense it has to be considered pri- 
marily as a working hypothesis which will be confirmed or 
contradicted in the future by more accurate data. 



2.3. Other experimental constraints 

A constraint of the amplitude of the fluctuation power 
spectrum at cluster scale can be derived from the cluster 
mass and the X-ray temperature functions. It is usually 
formulated as a constraint for the density fluctuation in 
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COBE2, Tegmark & Hamilton 1997 

COBE3, Tegmark & Hamilton 1997 

COBE4, Tegmark & Hamilton 1997 

COBE5, Tegmark & Hamilton 1997 

FIRS, Ganga et al. 1994 *' 

COBE6, Tegmark & Hamilton 1997 

COBE7, Tegmark & Hamilton 1997 

COBE8, Tegmark & Hamilton 1997 

Tenerife, Hancock et al. 1997 * } 

PythonVl, Coble et al. 1999 

iac/bartol2, Femenia et al. 1997 

SP91, Gundersen et al. 1995 

SP94, Gundersen et al. 1995 

Boomerang, Mauskopf et al. 1999 

BAM, Tucker et al. 1997 * } 

PythonV2, Coble et al. 1999 

QMap Fl+2Ka, de Oliveira-Costa et al. 1998 

MSAM, Wilson et al. 1999 

SKI, Netterfield et al. 1997 

Pythonl, Piatt et al. 1997 

Argol, de Bernardis et al. 1994 

Argo2, Masi et al. 1996 

Boomerang, Mauskopf et al. 1999 

PythonV3, Coble et al. 1999 

QMap F1+2Q, de Oliveira-Costa et al. 1998 

IAB, Piccirillio et al. 1993 *' 

QMap Fl+2Ka, de Oliveira-Costa et al. 1998 

TOC098, Miller et al. 1999 

MAX, Tanaka et al. 1996 *' 

Python V4, Coble et al. 1999 

TOC098, Miller et al. 1999 

Boomerang, Mauskopf et al. 1999 

SK2, Netterfield et al. 1997 

Python2, Piatt et al. 1997 

PythonV5, Coble et al. 1999 

MSAM, Wilson et al. 1999 

PythonV6, Coble et al. 1999 

Boomerang, Mauskopf et al. 1999 

TOC098, Miller et al. 1999 

PythonV7, Coble et al. 1999 

Sk3, Netterfield et al. 1997 

Boomerang, Mauskopf et al. 1999 

SK4, Netterfield et al. 1997 

Boomerang, Mauskopf et al. 1999 

TOC098, Miller et al. 1999 

SK5, Netterfield et al. 1997 

CAT1, Scott et al. 1996 * } 

Boomerang, Mauskopf et al. 1999 

MSAM, Wilson et al. 1999 

Ring5M2, Leitch et al. 1998 

CAT2, Scott et al. 1996 * ) 



*' - l e ff and band width were taken from Max Tegmark's CMB data analysis center (Tegmark 1999) 
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Fig. 2. Observational data of CMB fluctuation (Table 2) and a 
sixth order polynomial fit to a power spectrum (solid line) . The 
dotted lines restrict the space of fitting curves which deviate 
from best fit by less than la (A\ 2 = 47.9 for 44 degrees of 
freedom) . 
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Fig. 3. The contour of positions l p and amplitudes A p of first 
acoustic peak which corresponds to the range of fitting curves 
which are in the 68.3% range of probability of point distribu- 
tion. The box which contains ellipse gives la errors for l v and 
A p . The position l p and amplitude A p for best fit coefficients 
are shown as a cross (see also in Fig. 2). 



a top- hat sphere of 8h 1 Mpc radius, as, which can be 
calculated for a given initial power spectrum P(k): 



4 



i f°° 

= tt^ / fc 2 P(fc)W 2 (8M P c k/h)dk, 
2n Jo 



(3) 



where W(x) — 3(sinx — xcosx)/x z is the Fourier trans- 
form of a top-hat window function. A recent optical de- 
termination of the mass function of nearby galaxy clusters 
(Girardi et al. 1998) gives a 8 £l° m 46 ~°- 09n ™ = 0.60 ± 0.04. 
Several groups have found similar results using different 



methods and different data sets (for a comprehensive list 
of references see Borgani et al. 1999). To take into account 
the results from other authors we have decided to use more 
conservative error bars: 



cr 8 0' 



0.46-0.09Q„ 



0.60 ±0.08 



(4) 



^From the existence of three very massive clusters of 
galaxies observed so far at z > 0.5 a further constraint 
has been established by Bahcall & Fan 1998 



ct 8 0^ = 0.8 ±0.1 



(5) 



where a = 0.24 if Q A = and a = 0.29 if S1 A > with 
ft A + Q m = 1- The relation of this value to other tests will 
be analyzed too. 

Another constraint on the amplitude of the linear 
power spectrum of density fluctuations in our vicinity 
comes from the study of galaxy bulk flows in spheres of 
large enough radius around our position. Since these data 
may be influenced by the local super-cluster (cosmic vari- 
ance), we will use only the value of bulk motion - the 
mean peculiar velocity of galaxies in the sphere of radius 
bOh^Mpc given by Kolatt & Dekel 1997, 



V 50 = (375 ± 85)km/s. 



(6) 



An essential constraint on the linear power spectrum of 
matter clustering on small scales (k ~ (2 — 40)/i/Mpc 
comes from the Ly-a forest of absorption lines seen in 
quasar spectra (Gnedin 1998, Croft et al. 1998 and refer- 
ences therein). Assuming that the Ly-a forest is formed 
by discrete clouds with a physical size close to the Jeans 
scale in the reionized inter-galactic medium at z <~ 2 — 4, 
Gnedin 1998 has obtained a constraint on the value of the 
r.m.s. linear density fluctuations 



1.6 < a F (z = 3) < 2.6 (95%C.L.) 

at k F w 34f^/ 2 h/Mpc • 



(7) 



Taking into account the new data on quasar absorp- 
tion lines, the effective equation of state and the tem- 
perature of the inter-galactic medium at high redshift 
were re-estimated recently by Ricotti et al. 1999. As re- 
sult the value of Jeans scale at z = 3 has moved to 
k F w 380„ /2 VMpc (Gnedin 1999). 

The procedure of recovering the linear power spec- 
trum from the Ly-a forest has been elaborated by 
Croft et al. 1998. Analyzing the absorption lines in a sam- 
ple of 19 QSO spectra they have obtained the following 
constraint on the amplitude and slope of the linear power 

1/2 

spectrum at z = 2.5 and k p — 1.5flm h/Mpc, 



A 2 p (k p ) = fc3p(fc p )/27T 2 = 0.57 ± 0.26, 
. _ A log P(k) 



A log k 



-2.25 ±0.18, 



(8) 
(9) 
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S h = 1.95 x io _5 a 



-0.35-0. 19 In Q„ 



(13) 



(95% CL). In addition to the power spectrum measure- Universe without tensor mode and cosmological constant 
ments we will use the constraints on the value of the Hub- is given by 
ble constant 

h = 0.65 ±0.15 (10) 

For a flat model with cosmological constant (£l m +£l\ = 1) 
which is a compromise between measurements made by we have 
two groups: Tammann & Federspiel 1997 

and Madore et al. 1998. We also employ nucleosynthesis 
constraints on the baryon density of 



5 h = 1.94 X lQ-5^0.785-0.05 lnO me -0.95n-0.169ft 2 (u) 



n b h 2 = 0.019 ± 0.0024(95%CL) 



(n = n — 1). The normalization constant A is then given 
(11) by 



given by Buries et al. 1999. An earlier value of fl},h 2 = 
0.024 ±0.006 by Tytler et al. 1996 will be used to analyze 
the influence of this assumption on the obtained cosmo- 
logical parameters. 

3. Testing the Method 

In order to test our method to determine cosmological pa- 
rameters for stability, we have constructed a mock sample 
of observational data. We start with a set of cosmological 
parameters and determine for them the "observational" 
data which would be measured in case of faultless mea- 
surements with la errors comparable to the observational 
errors. We then insert random sets of starting parameters 
into the search program and try to find the right model 
which corresponds to the mock data. The method is sta- 
ble if we can recover our input cosmological model. Even 
starting very far away from the true values, our method 
reveals as very stable and finds the 'true' model whenever 
possible (see Table 3). 

One of the main ingredients for the solution for our 
search problem is a reasonably fast and accurate deter- 
mination of the transfer function which depends on the 
cosmological parameters. We use the accurate analytical 
approximations of the MDM transfer function T(k; z) de- 
pending on the parameters fl m , Q&, N v and h by 
(Eisenstein & Hu 1999 and Novosyadlyj et al. 1999). 

The linear power spectrum of matter density fluctua- 
tions is 



P(k; z) = Ak n T 2 (k; *)D?(z)/D?(0), 



(12) 



where A is the normalization constant and D\(z) is 
the linear growth factor, which can be approximated by 
(Carroll, Press & Turner 1992) 



Di(z) 



5 Q. m (z) 
2 1 + z 



1 209O TO (z) - Q* m (z) 4/7 ' 
70 + 140 + m { ) 



A = 2ttX(3000//i) 



3+n 



Mpc 4 



(15) 



The Abell-ACO power spectrum is related to the mat- 
ter power spectrum at z = 0, P(k; 0) by the cluster biasing 
parameter b c i. We assume scale-independent, linear bias: 



P A+ ACo(k) = b 2 cl P(k;0). 



(16) 



For a given set of parameters n, f2 m , f^, h, N v and 
b c i theoretical values of Pa+aco {kj ) can now be obtained 
for the values kj of Table 1. We denote them by yj ( j = 
1 13). 

The dependence of the position and amplitude of the 
first acoustic peak of the CMB power spectrum on cosmo- 
logical parameters has been investigated using CMBfast 
by Seljak & Zaldarriaga 1996. As expected, the results 
are, within sensible accuracy, independent of the hot dark 
matter contribution (f2„). This is illustrated in Figs. 4 and 
5. For the remaining parameters, n, h, fib and Ct\, we have 
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Fig. 4. The dependence of the acoustic peak amplitude A p on 
neutrino content f2„ 



where fl m (z) = n m (i + zf/ (fi m (i + z) 3 + fi A ). 

We normalize the spectra to the 4-year COBE data 

which determines the amplitude of density perturba- determined the resulting values i v and A v with CMBfast 

tion at the horizon crossing scale, 8h (Liddle et al. 1996, for a network of model parameters. The values £ p ,A p in- 

Bunn and White 1997), which for a matter dominated between grid points are then obtained by 4-dimensional 
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Fig. 5. The dependence of the acoustic peak position £ p on 
neutrino content Q. v 



interpolation. This allows a fast and sufficiently accu- 
rate calculation of the peak position and amplitude for 
a given set of parameters in the range 0.7 < n < 1.4, 
0.3 < h < 0.8, < n b < 0.2 and < fl A < 0.8 con- 
sidered in this work. The accuracy of this interpolation is 
estimated to be within 2%. We denote l v and A p by y\4 
and yi5 respectively. 

The theoretical values of the other experimental con- 
straints are obtained as follows: The density fluctuation as 
is calculated according to Eq. (3) with P(k; z) taken from 
Eq. (12). We set y w = a 8 n° m 4(i ' - 09n ^ and y 17 = a 8 fl a , 
where a = 0.24 for ft A = and a = 0.29 for fl A > 0, 
respectively. 

The r.m.s. peculiar velocity of galaxies in a sphere of 
radius R = 50/i _1 Mpc is given by 



1 r°° 

V 2 = — k 2 P (v) (k)e- k2R2 fW 2 (50Mpck/h)dk, (17) 
2tt Jo 

where pM(k) is power spectrum for the velocity field 
of the density- weighted matter (Eisenstein & Hu 1999), 
VF(50Mpc k/h) is the top- hat window function. A previ- 
ous smoothing of raw data with a Gaussian filter of radius 
Rf = 12/i _1 Mpc is employed here similar to the procedure 
which has led to the observational value. For the scales of 
interest pW(fc) w (n OM H Q ) 2 P(k;0)/k 2 . We denote the 
r.m.s. peculiar velocity by 

The value by Gnedin 1998 from the formation of Ly- 
a clouds constrains the r.m.s. linear density perturbation 

1/2 

at z = 3 and kp = 38£l„i ft/Mpc. In terms of the power 
spectrum of is given by 

a 2 F (z) = k 2 P(k; z^-^^dk, . (18) 

It will be denoted by yig. The corresponding value of the 
constraint by Croft et al. 1998 is 

A 2 p (k p ,z) = k 3 p P(k p ,z)/2ir 2 , (19) 



at z = 2.5 and k p = 0.008iJ(z)/(l + z)(km/s)-\ (where 

1 /2 

H{z) = H [Cl m (l + z f + n A ] is the Hubble parame- 
ter at redshift z) will be denoted by 7/20- The slope of the 
power spectrum at this scale and redshift, 



(20) 



A log P(k,z) 

n ( z ) = — t~i ; — i 

w A log k 

is denoted by j/21 • 

For all tests except Gnedin's Ly-a clouds we used 
the density weighted transfer function T c t> v {k,z) from 
Eisenstein & Hu 1999. For Gnedin's <r F we use T cb (k,z) 
according to the prescription of (Gnedin 1998). It must 
be noted that even in the model with maximal Q v (~ 0.2) 
the difference between T c b(k, z) and T c b„(k, z) is less than 
12% for k < k p . 

Finally, the values f2& and h are denoted by 2/22 and 
2/23 respectively. 

The relative quadratic deviations of the theoretical val- 
ues from their observational counterparts are given by \ 2 '- 



X 



23 

= E 

i=i 



Vi ~ Vj 



(21) 



where jjj and Afjj are the experimental data and their dis- 
persion, respectively. The set of parameters n, fi m , fib, h, 
Q, v , N v and b c i or some subset of them can be determined 
by minimizing \ 2 using the Levenberg-Marquardt method 
(Press et al. 1992). The derivatives of the predicted values 
with respect to the search parameters which are required 
by this method are calculated numerically using a relative 
step size of 10~ 5 with respect to the given parameter. 

The method was tested in the following way. Assum- 
ing a 4-year COBE normalized tilted AMDM model with 
the parameters n — 1.2, TO = 0.55, = 0.06, Vt v = 0.2, 
N v = 2, h = 0.65 and assuming further a cluster bias- 
ing parameter b c i = 3.0 we have calculated mock clus- 
ter power spectrum PA+ACo(kj) and treated them as ?7j, 
i = 1, ...,13. The remaining mock data jji, i = 14, ...,23 
have been calculated as described above. We have assigned 
to these mock data the same relative 'experimental' errors 
as in the corresponding experiments described in the pre- 
vious section. 

We then used these mock data to search the parame- 
ters n, fi m , fib, h, Q„, and b c i (N u was fixed). As starting 
parameters for the search program we assumed random 
values within the allowed range. We have searched for 
the parameters assuming the "true" value of two species 
of massive neutrinos as well as assuming three species of 
massive neutrinos. The parameters obtained for different 
cases are presented in Table 3. The errors in the deter- 
mined parameters are calculated as root square from di- 
agonal elements of the standard error covariance matrix. 
In all cases the code found all the previous known param- 
eters with high accuracy. This means that the code finds 
the global minimum of \ 2 independent of the initial values 
for the parameters. 



Cosmological parameters from LSS 



9 



Table 3. Test of the method: results of parameter search from mock data for the tilted AMDM model (n = 1.2, fl m = 0.55, 
fib = 0.06, fl„ — 0.20, N v — 2, h — 0.65). In test 1, all parameters are determined; in the 2nd to 6th tests, some parameters are 
fixed. For each test the first row corresponds to the case when number of species of massive neutrinos is equal the input value 
(2) and the second - when N v = 3. 



No 


N v 


y 2 ■ 

Amm 




n 






fl 




flv 






fib 


h 




b c i 




1 


2 


0.00 


1 


20±0 


07 





55±0.15 





200±0 


059 





060±0.022 


0.65±0.12 


3 


00±0 


45 




3 


0.05 


1 


21±0 


07 





63±0.17 





251±0 


073 





061±0.022 


0.65±0.12 


3 


12±0 


46 


2 


2 


1.72 


1 


18±0 


06 





79±0.07 





281±0 


055 





101±0.005 


0.50*) 


3 


58±0 


32 




3 


3.84 


1 


21±0 


06 





98±0.08 





446±0 


035 





101±0.005 


0.50*' 


3 


65±0 


30 


3 


2 


0.00 


1 


20±0 


07 





55±0.05 





200±0 


036 





060±0.003 


0.65*' 


3 


00±0 


27 




3 


0.05 


1 


21±0 


07 





62±0.05 





249±0 


043 





060±0.003 


0.65* 


3 


10±0 


27 


4 


2 


0.68 


1 


21±0 


07 





45±0.04 





168±0 


029 





045±0.002 


0.75*' 


2 


73±0 


24 




3 


0.80 


1 


22±0 


07 





51±0.05 





207±0 


034 





045±0.002 


0.75*' 


2 


83±0 


25 


5 


2 


0.00 


1 


20±0 


07 





55±0.05 





200±0 


036 




0.060*' 


0.65*' 


3 


00±0 


27 




3 


0.05 


1 


21±0 


07 





62±0.05 





249±0 


043 




0.060*) 


0.65*' 


3 


10±0 


27 


6 


2 


18.97 


1 


09±0 


06 




0.30" 





039±0 


003 





047±0.007 


0.73±0.05 


3 


60±0 


38 




3 


18.20 


1 


02±0 


08 




0.30*) 





000±0 


001 





064±0.021 


0.63±0.10 


3 


55±0 


31 



- fixed parameters. 



Our conclusions from the test results can be summa- 
rized as follows: 

1. If all parameters are free and N v = 2 (the input 
value) the code finds the correct values of the free param- 
eters (test 1, for N u = 2 in Table 3). 

2. If all parameters are free and N v = 3 the code finds 
values of the free parameters which are in the la range of 
errors (test 1, for N v = 3 in Table 3). 

3. If some parameters are fixed and differ from the 
input values (tests 2, 4, 6 in Table 3) the code finds for the 
remaining search parameters values close to the correct 
ones. The most stable and accurate value is fib- The results 
for n, £1„ and fl m are in the < 2a range of the correct 
values. The most uncertain solutions are found for n and 
fl v if an incorrect value for fl m has been assumed (test 6 
in Table 3). 

4. If some parameters are fixed to the predetermined 
ones and N v = 2 (the input value) the code finds the 
correct values of the free parameters (test 3 and 5 , for 
N v = 2 in Table 3), if N v = 3 the determined values are 
within the la range (test 3 and 5 , for N v = 3 in Table 3). 

In summary, the code determines the parameters n, 
fl u , fib, h, b c i and fl m correctly, if the observational data 
are correctly measured and the cosmological model as- 
sumed is correct; i.e. no curvature, a negligible amount of 
tensor perturbations and a primordial spectrum of scalar 
perturbations which is scale free from the present horizon 
size down to the scale of the Ly-a clouds. 



4. Results 

The determination of the parameters n, fl m , fib, h, fl v , 
N v and b c i by the Levenberg-Marquardt \ 2 minimization 
method can be realized in the following way: We vary the 
set of parameters n, flm, Qb, h, fl v and b c i or some subset 
of them and find the minimum of \ 2 . Since the N v process 
is discrete we repeat this procedure three times for N v =l, 
2, and 3. The lowest of the three minimums is the mini- 
mum of x 2 for the complete set of free parameters. The 
number of degrees of freedom Np = iV oxp — A?p ar = 7 if all 
parameters are free. It increases, if some of the parameters 
are fixed to a certain value. (Remember that even though 
we have 13 power spectra points, they can be described 
by just 3 degrees of freedom.) 

We have determined the minimum of \ 2 f° r N v =l, 2, 
3 in 11 different cases, where all observational data de- 
scribed in Sect. 2 are used. 

1) n, fl m , fl u , fib, h, and b c i are free parameters (Np — 

7); 

2) h = 0.5 is fixed, the remaining parameters are free 
(N F = 8); 

3) h = 0.6 (Saha et al. 1999,Tammann et al. 1999) is 
fixed, the remaining parameters are free (Np — 8); 

4) h = 0.72 (Madorc et al. 1998, Richtler et al. 1999) 
is fixed, the remaining parameters are free (Np — 8); 

5) h = 0.6 (Saha et al. 1999,Tammann et al. 1999) 
and h 2 fl b = 0.024 (Tytler et al. 1996) are fixed, the re- 
maining parameters are free (Np —9); 
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6) TO = 1.0 is fixed, the remaining parameters are free 
(N F = 8); 

7) Q TO = 0.3 is fixed, the remaining parameters are free 
(N F = 8); 

8) n = 1 is fixed, the remaining parameters are free 
(N F = 8); 

9) n = 1, TO = 1 are fixed, the remaining parameters 
are free (Nf — 9); 

10) n = 1, f2 TO = 0.3 and are fixed, the remaining 
parameters are free (N F = 9); 

11) Q„ = 7.4 x 10" 4 A„//i 2 is fixed by the lower limit 
of neutrino mass inferred by the observed neutrino oscil- 
lations in the Super-Kamiokandc experiment (Np = 8). 

For these 11 cases we find the minimum of \ 2 from 
which we determine the parameters presented in Ta- 
ble 4. Note, that for all models Xmin 1S m tnc range, 
N F - \/2Nf < Xmin < n f + V2Np which is expected for 
a Gaussian distribution of Np degrees of freedom. This 
means that the cosmological paradigm which has been as- 
sumed is in agreement with the data. (Note here, that 
the reduction of the 13 not independent data points of 
the cluster power spectrum to three parameters is very 
important for our analysis. Otherwise we would obtain a 
Xmin which is by far too small. If we would have assumed 
the 13 points of the Abell cluster power spectrum as in- 
dependent, resulting in Np = 17, the smallness of Xmin 
would have indicated that something is wrong in our ap- 
proach. But we might have drawn the wrong conclusion 
that the error bars be too large!) In Table 5 we present 
also the values of the different observational constraints 
for the best fit models found in Table 4. 




n=l & Q m =.3 fixed 



0.1 

k [h Mpo"'] 

Fig. 6. The observed Abell-ACO power spectrum (filled cir- 
cles) and the theoretical spectra predicted by tilted AMDM 
models with parameters taken from Table 4 (N v = 1). 

If all parameters are free (Table 4, case No 1), the 
model with one sort of massive neutrinos provides the 



best fit to the data, Xmin ~ 4-6- Note, however, that there 
are only marginal differences in X m in f° r ^ = 1? 2, 3. 
Therefore, with the given accuracy of the data we can- 
not conclude whether - if massive neutrinos are present 
at all - their number is one, two, or three. We summa- 
rize, that the considered observational data on LSS of the 
Universe can be explained by a flat AMDM inflationary 
model with a tilted spectrum of scalar perturbations and 
vanishing tensor contribution. The best fit parameters are: 
n = 1.12 ± 0.10, Q m = 0.41 ± 0.11, Q„ = 0.059 ± 0.028, 
N v = l,Q b = 0.039 ± 0.014 and h = 0.70 ± 0.12. The 
CDM density parameter is fl c dm = 0.31 ± 0.15 and I^a is 
considerable, ft a = 0.59 ± 0.11. 

The value of the Hubble constant is close to mea- 
surements by Madore et al. 1998. The spectral index co- 
incides with the COBE prediction. The neutrino matter 
density fi„ = 0.059± 0.028 corresponds to a neutrino mass 
m v = Mfl v h 2 w 2.7 ± 1.2 eV. The estimated cluster bias 
parameter b c i = 2.23 ± 0.33 fixes the amplitude of the 
Abell-ACO power spectrum (Fig. 6). All predictions of the 
measurements summarized in Table 5 are close to the ex- 
perimental values and within the error bars of the data. 

The predicted position of the acoustic peak (£ p = 215) 
is systematically lower than the experimental value deter- 
mined here from the complete data set on AT/T (£ p = 
253 ± 70). This position is nearly fixed by the requirement 
ft m + ft a = 1 and is only weakly dependent of the parame- 
ters varied in this study. The acoustic peak inferred by the 
Boomerang experiment (Mauskopf et al. 1999) is situated 
at £ <~ 200 and prefers models which are very close to flat 
(Melchiorri et al. 1999b). The models with low fl m ~ 0.3 
(case No. 7 in Table 4) fit the observable data somewhat 
less good than the best model (Ax min ~ 2.0) but all pre- 
dictions are still within the ler range. These models prefer 
a high Hubble parameter, h w 0.7 and no massive neu- 
trinos, fl v = 0. On the contrary, the matter dominated 
tilted MDM model (ft m = 1, models 6 in Table 4) prefers 
high Q„ = 0.22, three sort of massive neutrino and a low 
Hubble parameter, h — 0.47. This can be understood by 
considering one of the most serious problems of standard 
CDM, namely that the model, when normalized to COBE, 
has too much power on small scales. This problem can be 
solved either by introducing HDM and thereby damping 
the spectrum on small scales or by introducing a cosmo- 
logical constant which leads mainly to a 'shift of the power 
spectrum to the left'. 

Another interesting correlation can be seen in Table 
4, cases No 2-4, where we have fixed h. An increasing 
Hubble constant is compensated by a decreasing matter 
density, fi m , (i.e. increasing cosmological constant) and a 
decreasing baryon content due to the tight nucleosynthesis 
constraint on ftbh 2 . Furthermore, increasing the number 
of massive neutrino species N v from 1 to 3 leads to an 
increase of f2„ from 0.06 to 0.13 and to a decrease of f^A 
from 0.59 to 0.43 (case 1). 
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Table 4. Cosmological parameters determined for the tilted AMDM model with one, two and three species of massive neutrinos. 
In case No. 1 all parameters are free, in the other cases (No2-ll) some of them are fixed, as described above. 



No 


N v 


y 2 ■ 

Amm 




n 












Of, 


h 




b c i 




1 


1 


4.64 


1 


12±0.09 


0.41±0.11 





059±0 


028 





039±0.014 


0.70±0.12 


2 


23±0 


33 




2 


4.82 


1 


13±0.10 


0.49±0.13 





103±0 


042 





039±0.014 


0.70±0.13 


2 


33±0 


36 




3 


5.09 


1 


13±0.10 


0.56±0.14 





132±0 


053 





040±0.015 


0.69±0.13 


2 


45±0 


37 


2 


1 


7.50 


1 


11±0.09 


0.64±0.10 





075±0 


058 





076±0.005 


0.50*' 


2 


72±0 


28 




2 


7.46 


1 


12±0.09 


0.73±0.12 





120±0 


075 





076±0.005 


0.50*' 


2 


86±0 


28 




3 


7.46 


1 


13±0.09 


0.82±0.14 





163±0 


089 





076±0.005 


0.50*' 


2 


96±0 


29 


3 


1 


5.28 


1 


12±0.09 


0.51±0.07 





074±0 


041 





053±0.003 


0.60*' 


2 


43±0 


26 




2 


5.45 


1 


13±0.09 


0.59±0.08 





110±0 


053 





053±0.003 


0.60*' 


2 


56±0 


26 




3 


5.62 


1 


13±0.09 


0.66±0.10 





144±0 


063 





053±0.003 


0.60*' 


2 


66±0 


27 


4 


1 


4.67 


1 


12±0.10 


0.39±0.05 





058±0 


026 





037±0.002 


0.72*' 


2 


19±0 


23 




2 


4.84 


1 


13±0.06 


0.47±0.06 





101±0 


014 





037±0.002 


0.72*' 


2 


29±0 


18 




3 


5.12 


1 


14±0.10 


0.53±0.07 





130±0 


046 





037±0.002 


0.72*' 


2 


38±0 


25 


5 


1 


5.68 


1 


11±0.09 


0.53±0.07 





068±0 


043 




0.067*' 


0.60*' 


2 


49±0 


27 




2 


5.76 


1 


11±0.09 


0.61±0.09 





103±0 


056 




0.067*' 


0.60*' 


2 


62±0 


27 




3 


5.85 


1 


12±0.09 


0.67±0.10 





136±0 


067 




0.067*' 


0.60*' 


2 


71±0 


27 


6 


1 


12.23 


1 


07±0.09 


1.00*' 





116±0 


086 





118±0.027 


0.40±0.05 


3 


15±0 


39 




2 


10.18 


1 


10±0.09 


1.00*' 





177±0 


086 





099±0.022 


0.44±0.05 


3 


10±0 


38 




3 


8.80 


1 


12±0.09 


1.00*' 





219±0 


084 





085±0.019 


0.47±0.05 


3 


07±0 


38 


7 


1 


6.55 


1 


04±0.10 


0.30*' 





000±0 


005 





038±0.013 


0.71±0.12 


2 


25±0 


19 




2 


6.55 


1 


04±0.10 


0.30*' 





000±0 


005 





038±0.013 


0.71±0.12 


2 


25±0 


19 




3 


6.55 


1 


04±0.10 


0.30*' 





000±0 


005 





038±0.013 


0.71±0.12 


2 


25±0 


19 


8 


1 


6.21 




LOO* 


0.45±0.12 





042±0 


033 





038±0.014 


0.71±0.13 


2 


44±0 


31 




2 


6.60 




1.00*' 


0.50±0.14 





062±0 


043 





038±0.014 


0.71±0.13 


2 


57±0 


32 




3 


6.85 




LOO" 


0.51±0.14 





063±0 


012 





038±0.014 


0.71±0.13 


2 


69±0 


32 


9 


1 


13.01 




LOO*' 


1.00*' 





088±0 


075 





104±0.027 


0.43±0.05 


3 


23±0 


36 




2 


11.46 




LOO" 


1.00"° 





130±0 


073 





086±0.024 


0.47±0.06 


3 


23±0 


35 




3 


10.46 




1.00*' 


1.00*' 





159±0 


069 





075±0.021 


0.51±0.07 


3 


23±0 


35 


10 


1 


6.94 




1.00*' 


0.30*' 





000±0 


010 





034±0.009 


0.75±0.09 


2 


25±0 


20 




2 


6.94 




1.00*' 


0.30*' 





000±0 


010 





034±0.009 


0.75±0.10 


2 


25±0 


20 




3 


6.94 




1.00*' 


0.30*' 





000±0 


010 





034±0.009 


0.75±0.10 


2 


25±0 


20 


11 


1 


6.34 


1 


04±0.07 


0.35±0.10 




0.002* 


*) 





045±0.016 


0.65±0.11 


2 


36±0 


26 




2 


6.48 


1 


04±0.07 


0.36±0.10 




0.004* 


*) 





045±0.016 


0.65±0.11 


2 


38±0 


26 




3 


6.63 


1 


04±0.07 


0.37±0.10 




0.005* 


*) 





046±0.016 


0.64±0.11 


2 


41±0 


26 



*' - fixed parameters, **' mass density of neutrino is fixed by the lowest limit of the neutrino mass from Super-Kamiokandc 
results, h 2 Q.„ = \f&m?N„/MeV with VSm? = 0.07eV. 



If we use the nucleosynthe- 

sis constraint by Tytler et al. 1996 (case No 5), Xmin ^ s 
slightly higher than in case No 3. 

Now let us discuss models with a perfectly scale in- 
variant primordial power spectrum as predicted by the 
first inflationary models, n = 1 fixed (cases 8-10 in Ta- 
bles 4). If all of the remaining parameters are free (case 8) 
then this data set prefers a AMDM model with parameters 
n m = 0.45 ± 0.12 and h = 0.71 ± 0.13 and a somewhat 



lower neutrino content than the best fit model. Models 
with low matter content, fl m = 0.3, prefer a high Hubble 
parameter, h ~ 0.75 and no hot dark matter, f2„ = (case 
No 10). The matter dominated model TO = 1 (case No 
9) is the standard MDM model with Q u = 0.16 ± 0.07, 
three sort of massive neutrino (m, v = 1.3 ± 0.7eV) and 
h = 0.51 ± 0.07. 

If the HDM component is eliminated or fl u is fixed at 
the small value defined by the lower limit of the neutrino 
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Table 5. Theoretical predictions for the observational values of tilted AMDM models found with the parameters of Table 4 (for 
the value of N„ leading to the lowest x 2 ). 



No 


N v 


i v 


Ap 


_ oO-46-0.09n TO 




V5o,km/s 


<J F 


Ap(fc P ) 


Tip ( kp ) 


£o/10 9 yrs 


1 


1 


215 


84.1 


0.65 


0.74 


353 


2.03 


0.48 


-2.28 


12.3 


2 


1 


233 


91.6 


0.68 


0.71 


348 


1.87 


0.51 


-2.16 


15.0 


3 


1 


223 


87.8 


0.67 


0.73 


356 


1.91 


0.52 


-2.18 


13.5 


4 


1 


214 


83.2 


0.65 


0.74 


353 


2.04 


0.49 


-2.28 


12.2 


5 


1 


224 


89.3 


0.67 


0.72 


353 


1.89 


0.52 


-2.18 


13.3 


6 


3 


228 


84.7 


0.72 


0.72 


366 


1.78 


0.44 


-2.16 


13.9 


7 


1 


217 


80.7 


0.59 


0.70 


297 


2.08 


0.68 


-2.17 


13.4 


8 


1 


209 


69.5 


0.65 


0.72 


323 


1.88 


0.58 


-2.22 


11.9 


9 


3 


218 


68.9 


0.71 


0.71 


331 


1.72 


0.47 


-2.21 


12.8 


10 


1 


213 


73.4 


0.59 


0.70 


292 


2.00 


0.67 


-2.20 


12.5 


11 


1 


221 


80.9 


0.62 


0.71 


300 


2.03 


0.65 


-2.18 


14.1 


Obs. 


data 


253 ± 70 


80 ± 17 


0.60 ±0.08 


0.8 ±0.1 


375 ± 85 


2.0 ± .3 


0.57 ±0.26 


-2.25 ±0.2 


13.2 ±3 



mass J Sin 2 = 0.07 from the Super-Kamiokande experi- 
ment Q v = 7.4 x 10~ 4 Nv/h 2 , we obtain the best-fit value 
for the matter density parameter fi m w 0.39 ± 0.11 and 
Hubble constant h = 0.62 ± 0.12 (case No 11). 

The experimental Abell-ACO power spectrum and the 
theoretical predictions for some best fit models are shown 
in Fig. 6. Recently it was shown (Novosyadlyj 1999) that 
due to the large error bars, the position of the peak of 
P(k) at k w 0.05h/Mpc does not influence the determina- 
tion of the cosmological parameters significantly. Mainly 
the slope of the power spectrum on scales smaller than 
the scale of the peak position determines the cosmological 
parameters. 

The errors in the best fit parameters presented in Ta- 
ble 4 are the square roots of the diagonal elements of the 
covariance matrix. More informations about the accuracy 
of the determination of parameters and their sensitivity 
to the data used can be obtained from the contours of 
confidence levels presented in Fig. 7 for the tilted AMDM 
model with parameters from Table 4 (case No 1, N v = 1). 
The same contours for cases No 6 and 7 are shown in Fig. 8 
and 9, respectively. These contours show the confidence re- 
gions which contain 68.3% (solid line), 95.4% (dashed line) 
and 99.73% (dotted line) of the total probability distribu- 
tion in the two dimensional sections of the six-dimensional 
parameter space, if the probability distribution is Gaus- 
sian. Since the number of degrees of freedom is 7 they 
correspond to A% 2 =8.2, 14.3 and 21.8 respectively. The 
parameters not shown in a given diagram are set to their 
best-fit value. 



As one can see in Fig. 7a the iso-x 2 surface is rather 
prolate from the low-fi m - high-n corner to high-£! m - low- 
n. This indicates some degeneracy in n — £l m parameter 
plane, which can be expressed by the following equation 
which roughly describes the 'maximum likelihood ridge' 
in this plane within the la: 

nv/fU~= 0.73. (22) 

A similar degeneracy is observed in the fl„ — Q m plane in 
the range < VL V < 0.17, 0.25 < O m < 0.6 (Fig.7c). The 
equation for the 'maximum likelihood ridge' or 'degener- 
acy equation' has here the form: 

ft„ = 0.023 - 0.44O m ± 1.3f4 . (23) 

The 3rd column of Table 4 (Xmin) shows that all mod- 
els except 9th with N v = 1 are within the la contour of 
the best fit. 

The next important question is: which is the confi- 
dence limit of each parameter marginalized over the other 
ones. The straight forward answer is the integral of the 
likelihood function over the allowed range of all the other 
parameters. But for a 6-dimensional parameter space this 
is computationally time consuming. Therefore, we have 
estimated the la confidence limits for all parameters in 
the following way. By variation of all parameter we deter- 
mine the 6-dimensional x 2 surface which contains 68.3% 
of the total probability distribution. We then project the 
surface onto each axis of parameter space. Its shadow on 
the parameter axes gives us the la confidence limits on 
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cosmological parameters. For the best AMDM model with 
one sort of massive neutrinos the lcr confidence limits on 
parameters obtained in this way are presented in Table 6. 



Table 6. The best fit values of all the parameters with errors 
obtain by maximizing the (Gaussian) 68% confidence contours 
over all other parameters. 



parameter 


central value and errors 


n v 
n b 

h*) 

n 
b c i 


n 41+0-59 
u -^ 1 -0.22 

06+ 011 

o.o39i° :» 

70 +o.i5(+o.3i) 

19+0.27 
J - J " 4 — 0.30 

2- 22 lo.7 



* } - the upper limit is obtained by including the lower limit 
on the age of the Universe due to the age of oldest stars, to > 
13.2 ± 3.0 (Carretta et al. 1999). The value obtained without 
this constraint is given in parenthesis. 

It must be noted that the upper lcr edge for h is 
equal 1.08 when we marginalized over all other parameters 
and input observable data used here. But this contradicts 
the age of the oldest globular clusters t = 13.2 ±3.0 
(Carretta et al. 1999). Thus we have included this value 
into the marginalization procedure for the upper limit of 
h. We then have 8 degrees of freedom (24 data points) and 
the 6-dimensional x 2 surface which contains 68.3% of the 
probability is confined by the value 13.95. We did not use 
the age of oldest globular cluster for searching of best fit 
parameters in general case because it is only a lower limit 
for age of the Universe, besides it does not change their 
values as one can see from last column of Table 5. 

The errors given in Table 6 represent 68% likelihood, of 
course, only when the probability distribution is Gaussian. 
As one can see from Fig. 7 (all panels without degeneracy) 
the ellipticity of the likelihood contours in most of planes 
is close to what is expected from a Gaussian distribution. 
This indicates that our estimates of the confidence limits 
are reasonable. These errors define the range of each pa- 
rameter within which the best-fit values obtained for the 
remaining parameters lead to Xmin — 12.84. Of course, 
the best-fit values of the remaining parameters lay within 
their corresponding 68% likelihood given in the Table 6. 
It does however not mean that any set of parameters from 
these ranges satisfies the condition, Xmin — 12-84. 

For example, standard CDM model (Cl m = 1, h = 0.5, 
fib = 0.05, n = 1 and best-fit value of cluster biasing 
parameter b d = 2.17 (cr 8 = 1.2)) has Xmin = i42 (0: that 



excludes it at very high confidence level, > 99.999%. When 
we use the baryons density inferred from nucleosynthesis 
(h 2 fl b = 0.019 {b cl = 2.25, cr 8 = 1.14)) the situation does 
not improve much, Xmin = H2- Furthermore, even if we 
leave h as free parameter we still find Xmin = 16 (> 1c) 
with the best-fit values h = 0.37 and b c i — 3.28 (erg = 
0.74); this variant of CDM is ruled out again by direct 
measurements of the Hubble constant. 

The standard MDM model (fi m = 1, h = 0.5, fl b = 
0.5, ra = 1, fi„ = 0.2, N„ — 1 with a best value of the 
cluster biasing parameter b c i = 2.74 (cr 8 = 0.83)) does 
significantly better: it has Xmin = 23-1 (99% C.L.) which 
is out of the 2cr confidence contour but inside 3cr. With the 
nucleosynthesis constraint the situation does not change: 
Xmin — 22; also if we leave h as free parameter: Xmin — 21, 
h = 0.48. But if, in addition, we let vary fl v , we obtain 
x 2 min = 13 with best-fit values of fl v = 0.09, h = 0.43, 
b c i = 3.2 (cr 8 = 0.73). This means that the model is ruled 
out (as well as the model 9th in Table 4) by the data set 
considered in this work at ~ 70% confidence level only. 
But also here the best-fit value for h is very low. If we fix 
it at lower observational limit h — 0.5 then x m in = 18.9 
(the best fit values are: fl u = 0.15, b cl = 2.8 (er 8 = 0.83)), 
which corresponds to a confidence level of 95% . 

Therefore, we conclude that the observational data set 
used here rules out CDM models with h > 0.5, a scale 
invariant primordial power spectrum (n = 1) and Qk = 
n A = at very high confidence level, > 99.99%. MDM 
models with h > 0.5, n = 1 and Q k = Q A = arc ruled 
out at - 95% C.L. 

The best-fit parameters for 31 models which are in- 
side of lcr range of the best model arc presented in Table 
4. We conclude also that the observational data set used 
here does not rule out any of the 32 models presented in 
Table 4 at high confidence level but defines the lcr range 
of cosmological parameters for the AMDM models which 
match observations best. 

One important question is how each point of the data 
influences our result. To estimate this we have excluded 
some data points from the searching procedure. We have 
determined the best-fit parameters for the cases: 

— all points of Abell-ACO power spectrum Pa+aco (kj ) 
are excluded, 

— data on position and amplitude of acoustic peak, £ p , 
A p are excluded, 

— the 

value for cr 8 from Girardi et al. 1998, a s n m i6 ~ - om ™ 
is excluded, 

— the value for cr 8 from Bahcall & Fan 1998, cr 8 f2^ 29 is 
excluded, 

— both these tests are excluded, 

— the bulk motion, V50, is excluded, 

— the Ly-a constraint by Gnedin 1998 dp(z — 3) is ex- 
cluded, 
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b) 



d) 



0.80 r 



0.70 0.80 



0.70 0.80 



0.70 0.80 



0.70 0.80 



0.60 
0.55 - 

o.so : 

0.060 [ 



a 0.040 r 



0.060 f 



0.050 r 



a" 0.04 ; 



0.030 r 




1.30 r- 



0.10 

n 



0.15 0.20 




0.10 
Q 



0.15 0.20 




0.020 F .... i ... . 
0.50 0.55 0.1 



0.65 0.70 0.75 0.80 
h 




0.50 0.55 0.6 



0.65 0.70 0.75 0.80 
h 



Fig. 7. Likelihood contours (solid line - 68.3%, dashed - 95.4%, dotted - 99.73%) of the tilted AMDM model with N v = 1 and 
parameters from Table 4 (case 1) in the different planes of n — Q. m — Q v — ilt — h space. The parameters not shown in a given 
diagram are set to their best fit value. 
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Table 7. Parameters determined for the tilted AMDM with one sort of massive neutrinos if some of the data are excluded from 
the searching procedure. 



Excluded data 


\min 


n 




o„ 




h 


b c i 


All points of PA+Acoikj) 


1.45 


1.10±0.08 


0.42±0.11 


0.053±0.021 


0.041±0.013 


0.68±0.11 






4.22 


1.15±0.15 


0.40±0.11 


0.065±0.033 


0.040±0.016 


0.69±0.14 


2.20±0.33 


5 8 r^ 46 -°- 09n ™ 


3.68 


1.12±0.09 


0.47±0.15 


0.077±0.036 


0.042±0.016 


0.67±0.13 


2.19±0.33 




4.03 


1.11±0.10 


0.40±0.11 


0.052±0.031 


0.041±0.015 


0.68±0.12 


2.32±0.35 


Both cr 8 tests 


3.65 


1.12±0.10 


0.46±0.17 


0.072±0.048 


0.042±0.016 


0.67±0.13 


2.22±0.38 


Vso 


4.57 


1.11±0.10 


0.41±0.11 


0.057±0.030 


0.039±0.014 


0.70±0.12 


2.25±0.34 


of{z — 3) 


4.61 


1.13±0.11 


0.39±0.12 


0.056±0.030 


0.038±0.014 


0.71±0.13 


2.19±0.36 


A 2 p (k p , z = 2.5),n p (fep, z = 2.5) 


4.41 


1.13±0.10 


0.41±0.11 


0.069±0.035 


0.038±0.014 


0.70±0.13 


2.19±0.36 


Both Lya tests 


3.70 


1.11±0.10 


0.56±0.22 


0.222±0.291 


0.042±0.017 


0.67±0.13 


2.27±0.40 


h 


4.28 


1.11±0.10 


0.35±0.13 


0.051±0.024 


0.030±0.017 


0.79±0.22 


2.11±0.39 


n b h 2 


4.04 


1.14±0.09 


0.37±0.10 


0.068±0.017 


0.000±0.000 


0.66±0.10 


2.14±0.33 



— the Ly-a constraint by Croft et al. 1998 K 2 p {k pi z = 
2.5) and h p (k p , z = 2.5) are excluded, 

— both Ly-a tests are excluded, 

— data on the direct measurements of Hubble constant h 
is excluded, and 

— the nucleosynthesis constraint by Buries et al. 1999 is 
not used. 

The results for models with N v = 1 and all parameters 
free are presented in Table 7 (see for comparison model 
1 for N v = 1 in Table 4). Excluding any part of observ- 
able data results only in a change of the best-fit values 
of n, Ct m and h within the range of their corresponding 
standard errors. This indicates that the data are mutually 
in agreement, implying the same cosmological parameters 
(within the still considerable error bars). The small scale 
constraints, the Ly-a tests reduce the hot dark matter 
content from Q. v <~ 0.22 to <~ 0.075. The os-tests further 
reduce Q v to ~ 0.06. Including of the Abcll-ACO power 
spectrum in the search procedure, tends to enhance Q„ 
slightly. The most crucial test for the baryon content is 
of course the nucleosynthesis constraint. Its ~ 6% — In- 
accuracy safely keeps h 2 ^ near its median value 0.019. 
The parameter in turn is only known to ~ 36% ac- 
curacy due to the large errors of other experimental data 
used here, especially Hubble constant. The obtained ac- 
curacy of h (~ 17%) is better than the one assumed from 
direct measurements, ~ 23%. Summarizing, we conclude 
that all data points used here are important for searching 
the best-fit cosmological parameters. 

5. Discussion 

The best-fit parameters obtained in this paper are within 
the allowed range of parameters found by other authors 
using different constraints. For example, for the AMDM 
model with scale- invariant primordial power spectrum and 
one sort of massive neutrinos which contributes 10-20% of 



matter density, Valdarnini et al. 1998 found 0.45 < Q m < 
0.75 for h = 0.5, and 0.3 < ft m < 0.5 for h = 0.7. Similar 
constraints have been given by Primack & Gross 1998 for 
AMDM models with two species of massive neutrinos. Our 
values of Cl m are within these ranges. But at the boundary 
of this parameter range x 2 > 20, which is outside of the 
2a confidence contour. 



Recently Bahcall et al. 1999 have shown that the CMB 
anisotropy data, the cluster evolution and the SNIa 
magnitude-redshift relation indicate a flat Universe with 
accelerated expansion, compatible with a Oa — 0.7 and 
TO w 0.3 if CDM (Q„ = 0) is assumed. As we can see 
from Table 4 (case No 11, 1), our analysis leads to the 
same conclusion if we set density of hot dark matter to the 
minimum value compatible with the Super-Kamiokande 
experiment (less than 1% of the total density). However, 
if f2„ is a free parameter, the observational data consid- 
ered in this work lead to a AMDM with a slightly blue 
spectrum of primordial fluctuations (case No 1 in Table 
4). 

In our preferred tilted AMDM models (case No 1) the 
masses of neutrinos are vn v = 2.7 ± 1.2 cV for model with 
N v = 1, m v = 2.4 ± 1.0 eV when N v = 2 and m v = 
2.0 ± 0.8 eV for model with N v — 3. The accuracy of 
neutrino mass or density determination is modest because 
the observational constraints depend stronger on Sl m and 
n than on 0„ and N v . In models with fixed low matter 
density fl m = 0.3 (case No 7 and 10) the best-fit values 
of the neutrino density arc 0, i.e. even below the 

lower limit of the massive neutrino contribution to the 
cosmological density indicated by the Super-Kamiokande 
experiment. However, the la contours of the low Q models 
include the Super-Kamiokande limit (see Fig. 9b). 
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Fig. 9. Likelihood contours (solid line - 68.3%, dashed - 95.4%, 
dotted - 99.73%) of tilted AMDM with N v = 3, fixed VL m = 0.3 
Fig. 8. Likelihood contours (solid line - 68.3%, dashed - 95.4%, and parameters from Table 4 (case 7) in the differ ent planes of 
dotted - 99.73%) of tilted AMDM with N v = 3, fixed Q m = 1 n _ ^ _ ^ _ h gpace The parameters not shown in a given 
and parameters from Table 4 (case 6) in the different planes of diagram are set to their best fit value 
n — Q.„ — Q,h — h SDace. The parameters not shown in a eiven 
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In the last column of Table 5 we also indicate the age 
of the Universe, 



to 



3H 



1 



2n 



1/2 



(24) 



for each model as well as the age of the oldest globular 
clusters (Carretta et al. 1999). All models with parame- 
ters taken from Table 4 have ages which are in agreement 
with the oldest objects of our galaxy. 

We have used a scale independent, linear bias factor 
b c i as free parameter in order to fit the Abell-ACO power 
spectrum amplitude. 

Let us discuss in more detail how the model predictions 
presented in Table 5 match each observable constraint sep- 
arately. The predicted position of the acoustic peak for all 
models is lower than the one determined from the observa- 
tional data set presented in Table 2 (£ p = 253 ±70). Tilted 
AMDM models prefer l p ~ 210 - 230. This is due to the 
fact that the peak position depends only very weakly on 
the parameters discussed in this work. It is determined 
mainly by spatial curvature which we have set to zero 
here (together with the initial conditions which we have 
assumed to be adiabatic). However, our result is in good 
agreement with the most recent and so far most accu- 
rate determination of the peak position from one single 
experiment, the North American test flight of Boomerang 
(Mauskopf et al. 1999, Melchiorri et al. 1999b), which led 
to 0.85 < £l m + J1a < 1-25 with maximum likelihood near 
1 for adiabatic CDM models. The prediction of our best 
model for position of the first acoustic peak (£ = 215) 
matches the value given by Boomerang experiment I ~ 
200 very well. The central value from the combination of 
all available experiments, l p = 253, may very well be con- 
taminated by mutual calibration inconsistencies. 

Finally we want to discuss the possibility of using 
the averaged power spectrum of galaxies obtained by 
Einasto et al. 1999 to determine the parameters. This av- 
eraged spectrum of galaxies is determined in a wide range 
of scales (0.02/i/Mpc< k < 10/i/Mpc) and has substan- 
tially lower errors than the Abell-ACO power spectrum 
used here. Its la errors are <~ 4% on small scales and 
~ 20% at large scales versus <~ 40% and <~ 60% re- 
spectively for the Abell-ACO power spectrum. It is in- 
teresting to compare the predictions obtained from the 
power spectrum of galaxies with our analysis, because, 
as already mentioned in the introduction, the correc- 
tion of the linear power spectrum for nonlinear evolu- 
tion must be included into the algorithm. We use the 
fitting function by Smith et al. 1997, which transfers the 
linear into the nonlinear power spectrum, and the obser- 
vational constraint for the Hubble constant h — 0.6 ± 
0.02 (Saha et al. 1999,Tammann et al. 1999) as well as 
the nucleosynthesis constraint for the baryon content by 
Buries et al. 1999. Under these assumptions we find the 
following best fit parameters: n = 1.11 ± 0.02, h — 0.62 ± 



0.02, n m = 0.2±0.03, n b = 0.044±0.005, ft„ = 0.01±0.01, 
N u = 3 and galaxy biasing parameter b g = 1.52 ± 0.06. 

If wc add the remaining observations described in 
Sect.s 2.2 and 2.3 the best-fit parameters remain practi- 
cally unchanged due to the large number of (probably not 
independent) data points in the galaxy power spectrum. A 
model with these parameters has serious problems repro- 
ducing the experimental data set used here. Indeed, with 
these parameters we obtain x 2 m 61, for the data set used 
in the rest of this work, far outside 3a contour. The model 
predictions cr 8 f^ 46 -°- 09! ^ = 0.28 and a s Q° m 29 = 0.35 are 
<~ 4cr lower than the corresponding observational values 
by Girardi et al. 1998 and Bahcall & Fan 1998. Moreover, 
the peculiar velocity V50 is ~ 2a lower than the observed 
value, ap{z — 3) is and A^(k p ,z — 2.5) are ~ 3a and 
~ 1.5ct lower than the corresponding values inferred from 
the Ly-a measurements. Therefore, we conclude that a 
model with parameters determined by the galaxy power 
spectrum is ruled out by the observations discussed in this 
work. 

This result is not completely unexpected, because 
the galaxy power spectrum on small scales is probably 
influenced by a scale dependent bias (see for example 
Kravtsov & Klypin 1999, Fig. 3) which is not taken into 
account here. Moreover, the fitting formula for nonlinear 
evolution at k >1 h/Mpc may be incorrect. If we disre- 
gard the short wavelength part of galaxy power spectrum 
we find parameters close to those presented in Table 4. 



6. Conclusions 

Using Levenberg-Marquardt \ 2 minimization method we 
have determined the cosmological parameters of spatially 
flat, tilted AMDM models. We searched for a maximum 
of 6 parameters: the spectral index n, the matter content 
(^m + = 1), the hot dark matter content Q v , the 
baryon content Of,, the dimensionless Hubble constant h 
and the biasing parameter for rich clusters, b c i. The ex- 
perimental data set used in the search procedure included 
the Abell-ACO power spectrum (Retzlaff et al. 1998), the 
position and amplitude of the first acoustic peak in the 
angular power spectrum of CMB temperature fluctu- 
ations determined from the set of published measure- 
ments on different scales, the constraints for the den- 
sity fluctuation amplitude as derived from the mass func- 
tion of nearby and distant clusters (Girardi et al. 1998, 
Bahcall & Fan 1998), the mean peculiar velocity of galax- 
ies in a sphere of radius 50h _1 Mpc (Kolatt & Dekel 1997), 
the constraints on amplitude and tilt of the power spec- 
tra at small scales obtained from Ly-a clouds at z=2- 
3 (Gncdin 1998, Croft et al. 1998), the nucleosynthesis 
constraints (Tytler et al. 1996, Buries et al. 1999) and the 
COBE data (Bunn and White 1997) which is used to nor- 
malize the model power spectra. 

We have considered one, two and three species of mas- 
sive neutrinos. Wc have studied the influence of a reduc- 
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tion of the number of free parameters. In Table 4 we sum- 
marize the parameters which we have determined in 33 
different cases. Based on the results presented in Table 4 
we conclude: 

— The tilted AMDM model with one sort of massive neu- 
trinos and the best-fit parameters n = 1.12 ± 0.10, 
Q m = 0.41±0.11, 0\ = 0.59±0.11, Q. v = 0.059±0.028, 
tt h = 0.039 ± 0.014 and h = 0.70 ± 0.12 (standard 
errors) matches the observational data set best. The 
la (68.3%) confidence limits on each cosmological pa- 
rameter, obtained by marginalizing over the other pa- 
rameters, are 0.82 < n < 1.39, 0.19 < ft m < 1, 

o < o A < o.8i, o < n„ < 0.17, 0.021 < n b < o.rs 

and 0.38 < h < 0.85. 

— The degeneracies in the n — fi m and f2„ — f2 m planes 
n^/Tu = 0.73 and Q„ = 0.023 - 0.44f7 m + l.2,Vt 2 m are 
revealed. 

— For fixed Hubble constant h raising from 0.5 to 0.72, 
the best-fit value for il m decreases from 0.63 to 0.39 for 
AMDM models with N v = 1. For models with N v = 2 
and 3 the value of Vt m raises by ~ 0.08 and ~ 0.15 
respectively. The Q, v is higher for more species of mas- 
sive neutrinos, ~ 0.06 for one sort and ~ 0.13 for three, 
and decreases slowly for growing h. The inclination of 
initial power spectrum n correlates somewhat with £l v 
and grows slightly with h. 

— Fixing a low f2 m =0.3 a ACDM model without HDM 
matches the observational data set best. In this case 
the parameters are h = 0.71 ± 0.12, n = 1.04 ± 0.10 
and Q b = 0.038 ±0.013. 

— For all models the biasing parameter b c i of rich clusters 
is in the range 2.2-3.3, for the best model it equals 
2.23±0.33 (standard error). The la (68.3%) confidence 
interval is 1.5 < b c i < 3.5. 

— CDM models with h > 0.5, scale invariant primordial 
power spectrum n = 1 and J7a = = are ruled out 
at very high confidence level, > 99.99%. 

— Also pure MDM models are ruled out at <~ 95% C.L. 

Finally, we note that the accuracy of present observa- 
tional data on the large scale structure of the Universe 
is still too low to constrain the set of cosmological pa- 
rameters sufficiently, but we believe that our work shows 
the potential of such studies, which search for parameters 
including data from different, often complementary obser- 
vations. It is clear that with sufficiently accurate data, 
such a study may also reveal an inconsistency of model 
assumptions. 
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